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In this paper we revisit the glass model describing the 
macroscopic behavior of the High- Temperature superconduc- 
tors. We link the glass model at the microscopic level to the 
striped phase phenomenon, recently discussed widely. The 
size of the striped phase domains is consistent with earlier 
predictions of the glass model when it was introduced for 
High- Temperature Superconductivity in 1987. In an addi- 
tional step we use the Hubbard model to describe the micro- 
scopic mechanism for d-wave pairing within these finite size 
stripes. 

We discuss the implications for superconducting correla- 
tions of Hubbard model, which are much higher for stripes 
than for squares, for finite size scaling, and for the new view 
of the glass model picture. 

Keywords: High-T c , Superconductivity, Glass Model, 
Hubbard Model, Finite Size Scaling, Pseudo-Gap 



I. INTRODUCTION AND GLASS MODEL 

More then 10 years have passed since the introduction 
of the glass model in High- Temperature Superconduc- 
tivity |l] ||. The model was originally used for the de- 
scription of granular superconductors . But it became 
immediately clear that the glass behavior of the then 
new High-T c materials was not caused by the granular 
structure of the ceramic compounds; rather the magnetic 
fields used in the experiments suggested that the micro- 
scopic origin of glassy behavior occurred at a length scale 
smaller than that of the grains jl| . 

In this paper we revisit the glass model. Inspired 
by the recent study by Tsuei and Doderer JhJ of the 
charge confinement effect in cuprate superconductors, 
we link the intragrain inhomogeneities demanded in the 
lass theory to the microscopic concept of striped phase 
"ll,fl2| , and in particular to that of stripe domains [ p^|Jlo| 
which have recently become of wider interest to the High- 
Tc community. The stripe domains are identified with 
the areas of constant phase in the glass theory. Micro- 
scopically Tsuei and Doderer used their charge confine- 
ment model fiofl to explain various generic features of the 
normal-state pseodogap, including the magnitude of the 
gap. Such stripe domains may be described from first 
principles by a model such as the Hubbard Model, en- 
abling an ab initio numerical calculation of intra-domain 
pairing properties such as the effective pairing tempera- 
ture and pseudogap. 



In 0] a behavior reminiscent of spin glasses was de- 
tected, the quasi-de-Almeida-Thouless-line (quasi AT- 
line): 



H 2 ^ « T C (H) - T c (0) 



(1) 



where H is the magnetic field and T C (H) the correspond- 
ing critical temperature or to be more precise the tem- 
perature below which metastable behavior occurs as seen 
in [0] . The existence of the AT- line was first seen as an 
evidence for glass behavior, but it turned out, that a to- 
tally different physical mechanism was the cause for the 
H 2 / 3 behavior Jl]|. Numerical simulations of the High-T c 
glass model 0] lead to this conclusion. 

The glass model can be described by the Hamiltonian 

i 



H = - J C0S (<f>i - 4>j - > 



where 



<»j> 



Ay = I Adl 



A is the vector potential leading to Jl[] 



A i:i = —H 

00 



2-7T Xi + Xi 



(Vj - Vi) 



(2) 



(3) 
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In equation (g) to (^) </>o is the elementary flux quantum, 
H the magnetic field. <pi is the phase of the supercon- 
ducting wave function in an area or domain denoted by 
i; Xi and yi are the cartesian coordinates. Furthermore 
J is the coupling energy between two clusters according 
to |)) taken identical in the simulations in [0 , J? the line 
integral from cluster i to j, < i,j > the sum over the 
nearest neighbors, and finally Aij are the phase factors 
between cluster i and j. 

As a result of the simulations in (l) we consider fig- 
ure [l] which is reproduced from Jl| and follow-up papers. 
Because of the results of the simulations we had to con- 
clude, that the origin of the glass behavior lies inside the 
grains, therefore also inside single crystals. 

The if 2 / 3 behavior was only reproduced in two dimen- 
sions, therefore the domains or areas of constant phase 
4>i had to be located inside the CuO-planes, connected 
by weak links J (figure 0) . The size of the domains was 
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also known. Here we present a typical calculation of es- 
timating the domain size with data available at the time 
of the publication of [[l 14 1. 

As the magnetic field of reference we consider the up- 
per critical field if™. In figure || we see i?" in the phase 
diagram of equation (0), that H™ sw 0.5. We should de- 
note the 2tt/4>o was taken as unity as is normal in nu- 
merical simulations. 

For an estimation of the typical distance a between two 
clusters, we consider a plaquette in figure ||. 

In the simulation a = 1 . The idea is that the phase fac- 
tors are equal, Aij yS i m = A. 
H. (Ai jtSirn denoting A i3 
correspondingly in experiment) 
have Aij^ s i m = 0.5. Considering the experimental situa- 
tion we refer to figure 1 in (l4| . There the magnetization 
or magnetic moment fi(H ) is plotted versus H for various 
temperatures. Considering the lowest temperature n(H) 
tends to zero at 1.5T. Therefore H" w 1.5T (Tesla) 



ij,exp, for corresponding fields 
used in simulation and Aij^ exp 
Therefore for we 



in experiment with 



2.06 



lQ- 15 m 2 T we find that 
a 2 w 1 ■ 10 4 A^ and therefore a w lOOA. 

We wish to emphasize here, that the length scale of 
100 A deduced from the glass picture is consistent with 
the length of the stripe domains seen in STM data 
jjjl, and roughly with that of the stripe width deter- 
mined from neutron 11 and ultraviolet photoemission 
Ul data. We propose revision of the concept of clean, 
single crystal high temperature superconductor material 
to include the concept of striped phase-type inhomogene- 
ity with associated 'glassy' behavior. It should be noted, 
that the aging or memory effect in spin glasses was also 
detected 

At the end of this section we would like to counter a 
certain source of criticism concerning the glass model. 
Ten years ago the community only considered the case 
of s-wave superconductivity. In the model J = 1 was 
chosen. The glass behavior as described in fjj is totally 
dominated by frustration, which is caused by the .Ay's in 
equation (§). In s-wave superconductivity J would not 
change sign and therefore not lead to additional frustra- 
tion. This could be different in d-wave as pointed out by 
Sigrist and Rice jlBj . 

But the situation in the simulation of Q| was, that 
only weak " correlated" disorder was chosen in the frame- 
work of the square lattice. This means the coordinates Xi 
and yi in equation (^) were obtained by allowing a ran- 
dom displacement of the sites within a circle with radius 
r w 0.4 of the lattice constant around the original site 
of the square lattice. Therefore J does not change sign 
even in the case of d-wave superconductivity. It should 
be pointed out that only this type of weak "correlated" 
disorder (jfj in two dimensions reproduced the H 2 / 3 be- 
havior of the quasi AT-line jyj . At that time this result 
also lead to the prediction of a roughly two-dimensional 
plane-like anisotropy in the new superconductors |[| . 



Summarizing the glass model picture of fU is still ap- 
plicable to the correlated striped domains in the High-T c 
materials governing the d-wave symmetry and predicting 
the domain sizes correctly. 



II. MICROSCOPIC MECHANISM AND 
HUBBARD MODEL 

As pointed out the microscopic mechanism has to be 
considered inside the domains. While glass papers did 
not deal with the mechanism, it was clear at that time, 
that, if the glass picture applies, the observed T^ xp is 
T glass of the glasg mo del. The microscopic mechanism 
could easily have a much higher critical temperature T™. 

Inside the domains we propose an electronic mecha- 
nism as suggested by the Hubbard model (HM) or to be 
more precise the tt'-HM 



H — H-kin + 7~t 



pot 



with the kinetic 



E 



^i,j( C i,a- C j,o- c ja C iu) 



and the potential part 



n 



pot 



H,l 



(5) 



(6) 



(7) 



of the Hamiltonian. We denote the creation operator 
for an electron with spin a at site i with c\ , the corre- 
sponding annihilation operator with c i a , and the number 

operator at site % with n, iCT = cj a c i a . t{j is the hopping 
parameter only nonzero between nearest neighbors i,j 
(tjj = t) and next nearest sites i,j {t^j — t'), and finally 
U is the interaction being repulsive in the High-T c case. 

The sites of the HM can be associated with the copper 
in the CuO-layers of the High-T c crystals @ . 

We published evidence for superconducting behavior 
with d-wave symmetry of the repulsive HM as early as 
1994 Jl6|]. The d-wave symmetry in the superconduct- 
ing state was subsequently established for the High-T c 
materials fjj|. 

According to |2^^l| we use the (vertex) correlation 
function (instead of the largest eigenvalue) to provide ev- 
idence for superconductivity following the standard con- 
cept of off diagonal long range order |^|. The vertex 
correlation function (vertex CF) is defined for d x 2_ y 2- 
waves as 

Cj (r) =C d (r)-^J2 q (t, * + r)C [ (i + 6,t + r + 5') 

i,6,5' 

(8) 
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with the single particle correlation function C a = 
(cj a c i+r a ) for spin a to extract the pairing effects in 
the two particle CF 

Cd(r) = tJ2 9sgs'{c^4 +H c i+r+5tl c i+r1 ) (9) 

i,S,S' 

with the phase factors g s , g' s = ±1 to model the d x 2_ y 2- 
wave symmetry, the number of lattice points L and the 
sum 8 {8') over all nearest neighbors. 

In figure ^ we reproduce the evidence for d-wave pair- 
ing in the repulsive HM using the projector quantum 
Monte Carlo (PQMC) technique. These simulations 
[§HiH performed for the relatively weak-coupling value 
of U — 2, were found, with adequate numerical precau- 
tions, to be quite stable For higher [/-values, 
e.g. {7 = 4, PQMC is found to encounter difficulties 
for larger than 8x8 cluster sizes, and in fact evidence 
for finite plateaus was not reproduced, e.g. [§8|-§o|, for 
large U. However in the simulation technique of Zhang, 
Carlson and Gubernatis J3l| the existence of the plateau 
phenomenon in the vertex CF (figure was confirmed. 
Recent Variational Monte Carlo results |52| also confirm 
the existence of d-wave pairing in Hubbard clusters. 

Given the observation of plateaus supporting the ex- 
istence of d-wave pairing in finite clusters, what are the 
implications for bulk superconductivity? Up to now it 
was thought necessary to scale the finite cluster calcu- 
lations to infinite size in order to answer this question, 
which unfortunately presents difficulties as corrections to 
scalingare tremendously large compared to classical sys- 
tems |3^^,U] . However in view of the new stripe phase 
picture of the crystalline High-T c material, the simulated 
cluster can be reinterpreted as a single finite length stripe 
as seen in STM data |Q, with typical dimensions of 
roughly 12x4 unit cells, i.e. even smaller than our largest 
system sizes (16 x 16 in ]34]). Scaling to an infinite system 
is inappropriate as the PQMC cluster calculation is now 
descriptive of the actual quantum systems present in the 
crystal (although they are treated as totally isolated). 

We consequently performed additional simulations for 
12x4 systems and now report the remarkable result that 
the plateau (4 x 12) is about a factor of five higher than 
the plateau (12 x 12), figure ||. We averaged only the 
vertex CF CY(r) in the large range regime of r, i.e. for 
the distances \r\ > \r c \: 

c V = t- E c » ( 10 ) 

r, \T\>\r c \ 

with the number L c of points with \r\ > \r c \. The qualita- 
tive behavior of our results is not influenced by the choice 
of r c as long as we suppress the short range correlations 
(i.e. r c > 1.9). 

Due to finite size effects (see J|^,[54| or in more detail 
p6| ) the shape of the curves in respect to the filling in 
figure H is influenced by (n) and L. 



Considering the estimates of the microscopic criti- 
cal temperature T™ from grand canonical Monte Carlo 
(GQMC) |3ql simulations for the squared finite Hubbard 
model in fl33j we can calculate T™ from the x V,F = 
LC d ' versus /3 plot. (/? is the inverse temperature.) 
In figure ^ we reproduce the results of the GQMC runs 
in |33j. We take as estimate of T™ the inflection point 
of the Xd P versus (3 curve. From figure ^| we deduce 
(3 to be roughly (3 w 10. Taking the hopping param- 
eter t = 0.1 eV this corresponds to the temperature 
T c w 100K. 

In figure ^ all curves for different system sizes inter- 
sect at one temperature. Even so the magnitude of 
the superconducting signal at low temperatures is in- 
creased for striped geometries, we expect from figure 
that the value of T™ should be approximately the same. 
Here "T c " denotes the temperature where the correlation 
length equals the size of the systems. 

The first occurrence of preformed pairs in the finite 
Hubbard clusters is at an even higher temperature T* . 
T* is at the onset of the susceptibility, i.e. the tempera- 
ture with the first occurrence of a nonzero susceptibility, 
in figure [| A rough estimate is f3 « 5. With t = 0.1 eV 
this corresponds to T* « 200 K, which is in agreement 
with experimental results [|37H41[1. It should be noted, 
that the range for t is between 0.1 eV to 1 eV as deduced 
from the experiments l^- ^l, t hus this onset could be 
even at higher temperatures |45|. 

Recently Tsuei and Doderer |[o|], have pointed out a 
possible interpretation of the pseudogap as a function of 
doping in cuprates as the finite size gap in the quasipar- 
ticle spectrum of the stripe phase domains. Consistent 
with these ideas, we also identify the onset temperature 
T* with the pseudogap energy in the cuprate material. 
Numerical values for the energy gap between the ground 
state and the first excited state in finite size cluster sys- 
tems were obtained already in 1995 fj6j . Figure @ repro- 
duces a typical case for an y/8 x y/8 system as seen in 
Lanczos diagonalization (the finite size gap is large for 
the relatively small clusters for which Lanczos diagonal- 
ization is tractable). The experimentally seen supercon- 
ducting gap is different from the pseudogap (finite size 
gap). It is the gap seen in the macroscopic glass model 
when the system as a whole becomes superconducting 
and the finite gaps are correlated. Therefore it is clear 
that for the superconducting (glass model) gap a clear 
d-wave symmetry is seen [[L9| while in the case of the 
pseudogap the effect is smeared out (3^,0 as the exper- 
iment "averages" over different (finite size) clusters and 
slight correlated disorder. 
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III. SUMMARY AND CONCLUSIONS 

The new view of High-T c superconductivity emerging 
from Ref. Jl(]] and the present paper is based macroscop- 
ically on the interaction of nanoscopic, internally-paired, 
domains (glass model) and microscopically on the en- 
hancement of pairing found to occur within metallic do- 
mains constituting the striped phase. The glass model 
for example enables us to identify the superconducting 
T c with the Tf ass of the glass model. 

The microscopic model of a domain can be exploited, 
based on a Hamiltonian such as the Hubbard model, to 
estimate the relatively high temperature T* , the effective 
d-wave intra-domain pairing temperature, and to deter- 
mine the pseudogap which is identified [fiof with the finite 
size quasiparticle gap within the domain. In earlier sim- 
ulations we already found the first signs of pairing in the 
Hubbard model at remarkably high temperatures f§3| . 
We find in the present work a very large enhancement 
of the intra-domain pairing when a non-square cluster 
morphology, 12 x 4, appropriate for the striped phase, is 
considered. 

A further consequence of the glass model is that all 
changes in T£ xp are related to the interdomain coupling 
energy J of equation (^) and not to the microscopic mech- 
anism described by the Hubbard model. This means, 
e.g. for the isotope effect, that when masses are altered 
by changing O 16 to O 18 , leading to the isotope exponent 
a = AT C /AM, the mass effect is entering through an 
isotope shift in J or other glass model parameters. T* 
is probably not affected by the isotope effect ^7j. This 
would be another evidence for the combined glass model 

Hubbard model picture. 

In general we note that all experiments aiming at 
changing T™ and therefore the microscopic mechanism 
have to be reinterpreted or redone in terms of J of the 
glass model. Future research could usefully concentrate 
on these weak links as seems to be occurring fig}] . Only 
the glass model J governs the T% xp (in fact T c m J), 
and additionally the critical currents of the materials. 
Progress in increasing T c or the critical current so far 
has been only accidental as the underlying glass behav- 
ior was not recognized. 

In the scientific discussion about ten years ago it was 
concluded @ that the glass model would be only appli- 
cable to so called "bad" or "glassy" samples, i.e. samples 
with weak links described by Hamiltonian of equation 
& 01 > At this time future research had to concentrate 
on "good" samples without weak links. Here we argue 
again, that there are no "good" samples, as the T^ xp, s 
for all samples, single crystals, thin films, and ceramics, 
are for a definite composition about the same. It fol- 
lows consequently that the T£ xp is the Tf ass of the glass 
model. 

Furthermore it is clear as we deal with finite size do- 



mains (or stripes) that mean field theory and therefore 
BCS does not apply. 

Future research could usefully emphasize the origin 
and manipulation of the weak links (i. e. the J in the 
glass model) between the domains (or stripes). 
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FIG. 1. Schematic illustration of the physical situation in 
a single crystal of a High-T c ceramic, (analogous to figure 1 
in [7]) 
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FIG. 2. Phase diagram in H — T plane. Reversible 
weak-random (WR) phase between H ~ 0.0 and H l c ?s 0.01 
- 0.03. Between H c and w 0.5, superconducting glass 
phase. Irreversible effects are separated from reversible ef- 
fects by dashed line. Above H™, usual XY-spin glass phase 
(SG). Dashed line indicates freezing transition, (analogous to 
figure 10 in [1]) 
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FIG. 4. Vertex-CF versus Cooper pair distance \r\ for 
extended s-wave (xs) and d x 2_ y 2-wave (d) for a 8 x 8 x 1 
system with U = 2, t' = —0.22 and (n) = 0.78. (Inset shows 
the long range regime enlarged.) 



As 



2 
1 





-1 
0.65 



x 10" 3 






12x4: g g 






12x12: x- - -x 




x- 


x- -x 



0.7 0.75 0.8 0.85 0.9 0.95 



<n> 



Frustration 

FIG. 3. Plaquette with domains of constant phase (pi and 
phase differences Ay (see text) 



FIG. 5. Averaged (jr c | = 1.9) vertex CF versus filling 
for the tt'-HM with U = 2, t' = —0.22 and the system sizes 
12 x 4, and 12 x 12. The parameter of the PQMC were = 8 
and t — 0.125. Lines are to guide the eye. 
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FIG. 6. Cumulated d-wave Vertex-CF x d ' versus inverse 
temperature (3 for the two-dimensional Hubbard model with 
U = 2, t' = -0.22 and (n) « 0.8. (data from figure 5.10 in 
[33].) 
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energy 

FIG. 7. Density of states (DOS) for the two-dimensional 
Hubbard model with U = 2, t' = -0.22 and (n) w 0.75 and 
L = 8. Energy axis in units of t. (data from [46]) 
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